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Quasi-resonant scattering of light in two dimensions can be described either as a scalar or as a 
vectorial electromagnetic wave. Performing a scaling analysis we observe in both cases long lived 
modes, yet only the scalar case exhibits Anderson localized modes together with extremely long 
mode lifetimes. We show that the localization length of these modes is influenced only by their 
position, and not their lifetime. Investigating the reasons for the absence of localization, it appears 
that both the coupling of several polarizations and the presence of near-held terms are able to 
prevent long lifetimes and Anderson localization. 


I. INTRODUCTION 

Multiple scattering of waves has been the subject of 
intense debates in the context of disorder-induced An¬ 
derson localization [T]. Indeed, since the proposal to use 
electromagnetic waves in random media instead of elec¬ 
trons in solids exploiting the non-interacting properties of 
photons at low intensities [5] , many experiments and the¬ 
oretical studies have been performed. However, despite 
a decade-long research, the mere existence of Anderson 
localization of light [Si] and its relation to another long 
predicted phenomenon, namely Dicke super- and subra¬ 
diance |5] , are still not clearly understood |7| . The advent 
of laser-cooled atoms and their use to study both localiza¬ 
tion and super- and subradiance motivated the develop¬ 
ment of ab initio models of interference effects in multiple 
scattering of light |S] . As most experiments are typically 
performed in a three-dimensional setting, models have 
also been focused on such 3D configurations. However, 
both numerical and fundamental aspects of localization 
strongly depend on the dimension of the explored sys¬ 
tem |9]. For this reason we have focused our efforts on a 
2D system, where a precise study of the eigenvalues and 
eigenmodes of the system is more efficient than in 3D, 
because larger ’volumes’ can be simulated for a given 
number of scatterers. One further advantage is that the 
reduced dimensionality allows for a direct comparison of 
eigenvalues and eigenvectors between two regimes of scat¬ 
tering, one of a scalar model of light, the other of a vec¬ 
torial model of light where the wave polarization needs 
to be accounted for. This comparison recently revealed 
important differences observed for the eigenvalues of the 
relevant effective Hamiltonians HEiiin]. 

In this work, we investigate resonant scattering in 
a two-dimensional set-up, i.e., the light scattering and 
propagation are confined to two dimensions. This config¬ 
uration may be realized, e.g., with a disordered arrange¬ 
ment of scatterers in microwave cavities m , in photonic 
crystals near surface plasmons or with laser-cooled 
atoms located in an off-resonance optical cavity. In this 
geometry, the polarization orthogonal to the plane, called 
s-polarization, cannot couple through the scatterers to 


the planar (or p-) polarizations, hence it is described 
by a scalar light model. The two p-polarizations of the 
electromagnetic waves do couple, and this vectorial-like 
scattering includes near-held terms (see scheme in Fig0. 
Rotating the polarization of an incident wave allows to 
switch between the scalar or vectorial regime, between 
the presence or the absence of polarization degrees of 
freedom and near held terms, making it an ideal tool 
to investigate the role of polarization in localization and 
subradiance. 

In section H we present a detailed derivation of the lin¬ 
ear differential equations that rules the population evo¬ 
lution of atomic transitions. The spectral properties of 
these equations are investigated in section HI, where 
scalar and vector scattering are compared via scaling 
analysis. Finally, we present our conclusions in section IV 
highlighting uncorrelation between spatial and temporal 
localization of light. 



FIG. 1. (Color online) Two-dimensional scattering scheme: 
The radiation, of wavenumber k close to the atomic transi¬ 
tion ka = LUajc K, k, is Confined in the {x,y) plane, i.e., it has 
wavevectors of the form k = A:(cos 0, sin 0, 0). Two eigenvec¬ 
tors are shown: a localized s—polarized mode in the upper 
right part and a extended p—polarized mode in the left part. 
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II. MODEL 

Two-dimensional light scattering was investigated in 
microwave cavities, where light with a polarization or¬ 
thogonal to the plates approximately obeys Helmholtz 
2D scalar equation, and Anderson localization was ob¬ 
served m- Another possibility to emulate 2D light scat¬ 
tering is a cloud of cold atoms with no Doppler broaden¬ 
ing located inside an off-resonant optical cavity made by 
two metallic disks whose diameter is much larger than 
their mutual distance. This system, which is closer to 
situations previously studied in 3D, constitutes the toy 
model system we will explore in this paper. Let us 
consider an homogeneous disk-shaped cloud of N mo¬ 
tionless atoms sitting at randomly distributed positions 
Tj = (xj,yj,Zj) with j = I,-- - ,N, for which non- 
radiative interactions are neglected. Instead, only virtual 
and real photons couple the atoms within the optical cav¬ 
ity (axis z) whose resonance frequency is signihcantly de¬ 
tuned from the atomic transition uja- The electric dipole 
transitions occur between one non-degenerate ground 
state Igj), related to angular momentum £ = 0 and a 
triply degenerate excited state |e™), where m = 0, ±1 
indicate the projections of the angular momentum £ = 1 
over the quantization axis z. We consider 2D scattering 
restricted to a radial direction in the {x, y) plane with a 
surface density of the atomic cloud p = N/ttR'^, where R 
is the cloud radius. 

The interaction of the atoms with the radiation field is 
given by the following Hamiltonian: 

j=l k.s ^ 

j — l m— — l k,s 

with u!a the atomic transition frequency, CTzj the atomic 


diagonal term, whereas = \gj) and = 

{gj\ are the lowering and lifting atomic operators. 
dj(. ^ and dk.s refer to the creation and annihilation of a 
photon for mode k, with frequency Wfc . The coupling co¬ 
efficient reads g^'^ = ek,s ■ dj_m\/wk/2deoH, with V the 
quantization volume and = {gj \ erj \ej^ra) the dipole 
matrix element (with e is the electron charge). Thus, 
the two first terms in Q are the free energy contribution 
and the last term corresponds to the interaction with the 
vacuum modes. 

We then use the commutation relations 
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where the approximation in the first of these equations 
correspond to the linear optics regime |ej_m) — 

\gj) {gj\ « ("''^6 eliminate the possibility of multi¬ 

excitation in the system), to obtain the Heisenberg equa¬ 
tions for the operators: 
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dt 
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Equations ( |3a[ ) and ( |3b| ) show a correlated dynamics be¬ 
tween atomic levels and vacuum modes. 


The radiation field plays the role of a reservoir for 
atoms and is composed of an infinite number of degrees 
of freedom k, s, so it is convenient to trace over these. 
Using the unitary transformations dj™^ —>■ and 

dk,s —> dk,se*“’‘‘, we obtain the reduced equation evolu¬ 
tion for the atomic open system 


dt 


\ ^ \ ^ f J ( l,n '-(n) f, \ iuj„T , I, 

EE^^k’« +9k. 






l.n k.s 


( 4 ) 


We can apply the rotating wave approximation and ne¬ 
glect the fast oscillating terms proportional to As¬ 

suming that the photon transit time inside the atomic 
cloud is much shorter than the emission decay time, we 
can perform the Markov approximation dj"^ (t — t) « 
d;"^ (t) so the atomic transitions evolves according to the 


closed set of equations 


d^(rn) ^ 1 N 
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where the scattering kernel is defined as 
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with Tji = Vj — r;, and where the upper limit in the above 
integral has been extrapolated to t —)■ oo according to the 
Markov approximation. The spontaneous emission pro¬ 
cesses is naturally three-dimensional and the above scat¬ 
tering kernel a priori contains all light modes in 3D space, 
so all transitions may be coupled. Under the assumption 
of an effective two-dimensional scattering of light, we can 
perform in the continuous 3D density of modes the fol¬ 
lowing approximation: 

rf^k5(0-7r/2), (7) 

where 6 corresponds to the azimuthal angle in spheri¬ 
cal coordinates, and V now refers to the quantization 
volume delimited by the cavity. Observe we are treating 


the radiation field inside a volume, however the density of 
modes is practically parametrized in polar coordinates as 
k = fc (cos (j), sin </), 0). The usual relation for polarization 
vectors J2s ®k,s®k.s = ~ k/ik,^ (/r, i' the Cartesian 

components) here turns into 

+ E - k^k, 

Therefore, the 2D condition Q decouples the component 
of the dipole matrix elements from and dJ„j. 

This phenomena is absent in 3D scattering, where all 
components of dj^m are coupled, whereas at the other 
end, ID case exhibits all components trivially uncoupled. 

The single atom decay into two-dimensional vacuum 
modes is given by 


E 7 ,m* Ln 

9i,s 5k,s = 


kc 


2heoV 


A (r,7 = 0) = d-^df,. 


ck^dk 


lo {2Tr) heo Jo 


d<t) / ( 

Jo 
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where the index 2D applied on means the general 
kernel particularized to the two-dimensional scattering. 
The time integral present in equation solves by using 
the relation 


roo 

/ =-(5(/c±fca)±iP— 

Jo C UJa 


1 


± OJk 


( 10 ) 


where P refers to the Cauchy principal value. This term 
gives rise to the Lamb shift, a single atom energy shift 
which is due to its interaction with the radiation field. 
We will here neglect it, as it simply corresponds to a 
renormalization of the energy; remark that we do not 
neglect the so-called collective Lamb shift, that rises from 
the interaction between the atoms via virtual photons, as 
it is still present in the final scattering kernel. 

Using the relations 

= , (11a) 

I ^^m,n k^kj^^ = yi, (11b) 


where the expression for dj „i in (Hb) includes choosing 


the quantization axis over z, and plugging (101 into (§, 
we obtain the two different decay rates 


To — d^o,o {^ji — 0) — 


k^d'^ 


Ti = K±i^±i {vji = 0) = 


Airheo 
kU^ 


(12a) 


T, 


Sirheo 


= —. (12b) 


By looking at (10), one can see only the term k = ka 
contributes. Therefore, the uncoupling of d| ^ from the 
other components of dj ^ causes an anisotropy on spon¬ 
taneous emission process, since the lifetime of transition 
Cj,m —0 t Qj is twice shorter than —>■ gj. These 

decay rates will predict the coexistence of two scattering 
subsystems with different time scales, which vector na¬ 
ture of light will be crucial to select each subsystem is 
active. 

Finally, we address the collective term by calculating 
the integrals describing the coupling between the atoms 
via the radiation field 


A (r,7 0) = d-^df^„ / , 


r del, (e—_ c.c) 
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(13) 


Since the density of modes is non-zero only in the (a;, y) plane, the z component of the atoms positions does not 
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come into play, so in the relation 


k in the integral as ka- Using the relation 


k^k,^e 


,± 2 k-r. 


ji —_ 


dx^;,dx-; 




(14) 


/ dkJo {krji) sin (fccr) = —Hq (karjH^) 
Jo 


where Xji actually spans only (xji,yji). The angular in¬ 
tegral then reads 



^±ikx,t cos 4>±tky„ sin0^^ ^ 2tt Jq (krji) , 


(15) 


where Ha is the Hankel function of the first kind and of 
order a, we can calculate {xji ^ 0) from the action 

of the second order derivative of Hq (kaTji) with respect 
to Xji- Practically, we get 


where Jq denotes the Bessel function of the first kind and 
of order 0, and Xji = 


4 + 


the Euclidean distance 

between each pair of atoms in the plane. Despite the in¬ 
tegrands in (13) diverge in the limit fc —>■ oo, the modu¬ 
lus of the wavevectors k vary only slightly around k = ka 
(quasi-elastic scattering). We then apply the Wigner and 
Weisskopf approximation which approximates powers of 


^ 0) = Poido (karji ), (17a) 

^±p±i ^ 0) = Piifo ikaTji ), (17b) 

^ 0) = ^ 1-^2 (karji) . (17c) 

These coefficients allow to obtain the following set of 
equations for the atomic operators (t4^: 


dt 


dt 


2 ' 


ho 
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Z=1 



(± 1 ) 

j 


Y (-^0 {krji) , 

1^1 


(18a) 

(18b) 


wher e tan i pji = {yj — yi)/{xj — xi). In both equations 
(18a I and (18b), the atoms are coupled together through 
the same sublevel with a kernel term Ho^kr) that scales 
as a 2D spherical wave j^/r at large distances. This 
l/Y scaling in 2D corresponds to long range coupling 
so we expect global coupling (or cooperative effects) to 
dominate over nearest neighbor coupling. In the vecto¬ 
rial case, the m = ±1 sublevels are additionally cou¬ 
pled via a H 2 term which also scales as at long 

range. However, with respect to global versus local in¬ 
teractions, Hq diverges at the origin only as log(/cr), so 
that the contribution of a small volume around the par¬ 
ticle is finite {Jq° log{kr)2Trrdr < 00 ); Instead, the H 2 
term diverges at the origin as I/r^, which means that 
the interaction between the ±1 sublevels is dominated 
by the close neighbors/near-field terms at high densities 


(J^ ° 2'Krdr/r^ 


r_—>-0 


—27r log(r_)). 


Differently from the 3D case where scalar light is only 
an approximation for dilute systems, and where all sub- 
levels are normally coupled, in 2D geometries scalar 
model holds for high densities. Yet, controlling the polar¬ 
ization of the injected light allows to select either purely 
scalar or vectorial properties which make our approach 
quite versatile. In the end, two decoupled scattering sub¬ 
systems appear: one involving a single sublevel of the 
excited state (the scalar case), the other one involving 
the remaining two sublevels (the vectorial case). As one 
can note the scalar and vector kernels are not decoupled 


through energy shifts like 3D work in Ref. [5], but rather 
by geometrical constraints. The microwave or optical 
cavity reshapes the density of electromagnetic modes into 
two dimensions. 

The scattered field at a point r = (x, y) is calculated 
by a superposition of annihilation operators, namely 

E (r) = ^ ekek,sak..e-*--‘*+*‘^(19) 

k,s 


where £k = \/IiWk/2eoU. With similar procedures used 
up to here, using (3b), it can be shown to lead to the 
following classical equation 


ST, 
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2di 

HTi 


/ = 1 
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® I"" “ u|) + 


HO) 


T^ka |r - nl 

5] ^ e™e-2™^'i72(fcjr-r,|)d4. 


I—I m—dzl 
N 


X emHoika\r-ri\)al 


( 20 ) 


m=±l 


with taaipi = {y — yi)/{x — xi) and e± = (x± iy) /\/2. 
The m = 0 sublevel is coupled only to light with polariza¬ 
tion along z, whereas the m = ±1 sublevels are coupled 
together through p-polarizations. 
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III. SPECTRAL ANALYSIS OF THE LINEAR 
EQUATIONS 


We now turn our attention to the spectral properties 
of the system. The scatterin g mo des a re th e eigenmodes 
of the linear equations (18a) and (18b), where n la¬ 
bels the eigenmodes. Their lifetime l/7n and energy 
are given by the real and imaginary part of the associ¬ 
ated eigenvalues, respectively. Defining 
the modes can also be characterized by their inverse par¬ 
ticipation ratio (IPR) that 


quantifies the (inverse) number of atoms substantially 
involved to the scattering mode. In the vectorial model 
we renormalize the IPR to remain 1/2 for pairs. 



FIG. 2. (Color online) Inverse participation ratio of the scat¬ 
tering modes in the complex plane of the eigenvalues (7n,a;n) 
for scalar light (a) without and (c) with near-field terms, in 
units of Fo, and vectorial light (b) without and (d) with near- 
held terms, in units of Fi, (a) and (d) being the physical 
cases. Simulations realized for an homogeneous disk cloud of 
N = 5000 particles with an homogeneous p/k^ = 1 density. 


In the scalar case, the eigenvalue distribution shown 
in Figj^a) exhibits strongly subradiant modes, which we 
define as modes with very long lifetimes ( 7 „ ^ Pq). The 
distribution can be used to look for a single parameter 
scaling, by computing a spectral overlap function con¬ 
veniently defined as g = /{<-an — w„_i), where 

the modes n are ordered by increasing energy. In line 
with the 3D results, we observe a monotonic decrease 
of g with the system size for scalar light (see Fig|^a)). 
Consequently, the scaling function /3 = d\ag/ d\^kR) 
is clearly negative for all values of g (see Fig|^b)), as 
expected for Anderson localization in 2D. We note that 
this function g is only one among several possibilities of 
defining a spectral overlap and has not been shown to be 


unequivocally related to transport properties of electro¬ 
magnetic radiation. 

The dimensionless scaling parameter g is defined in the 
scaling theory m as the ratio between the Heisenberg 
time and the Thouless time. The former corresponds to 
the time associated to the mean spacing between the en¬ 
ergy levels, i.e. h/{En — i?„_i), which in our case reads 
1 /{ujn — U!n-i)- The latter corresponds to the time nec¬ 
essary to a photon to escape from the sample, and in our 
open system with eigenmodes of lifetime 1/jn, we de¬ 
fine it as (l/ 7 „). Following this interpretation, the local¬ 
ization regime is characterized by a Thouless (diffusion) 
time that becomes larger than the Heisenberg one. 



FIG. 3. (Color online) (a) Thouless number ^ as a function 
of the dimensionless system size kR, for density p/k^ = 0.3 
(b) scaling function /?• Longest mode lifetimes for different 
number of particles (c) as a function of the optical thickness 
bo and for all four scattering models and (d) as a function of 
the dimensionless parameter N^^^/kR for scalar light. The 
circles are for scalar light, the squares for vectorial light, the 
stars for scalar light with near field and the plus for vectorial 
light without near field. The black line in (c) stands for the 
~ l/6o curve of the radiation trapping lifetime. In (c) and 
(d), blue, brown and magenta curves correspond respectively 
to simulations with N = 1000, 2000 and 4000 particles. The 
insert in (d) shows Tmin as a function of atomic density, indi¬ 
cating that density is not the best scaling function for F^m, 
even in the high density limit. 

In the vectorial case exhibited in Figj^d) we observe 
a dramatically different eigenvalue distribution. Indeed, 
even though long-lived modes ( 7 „ <C Fq) exist, they are 
limited to values larger by several orders of magnitude 
compared to the scalar case. The corresponding spectral 
overlap also shows a distinct behavior, with g almost in¬ 
dependent of the system size kR [see Figj^a)], yielding 
a scaling function (3 close to zero, albeit slightly negative 
[see Figj^b)]. This behavior of the scaling function fi 
might make this vectorial case very interesting to study 
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fine corrections of the atom-atom interactions as it seems 
to be close to the critical regime. 

The above discussion is consistent with the conclusions 
drawn from the study of eigenvalues in 3D Ulin]. With 
the aim to pin down the essential ingredient of the dif¬ 
ference between the scalar and vectorial model, we arti¬ 
ficially introduced or removed short range terms in the 
two configurations. More specifically, we removed the 
near-held coupling from vectorial scattering by substitut¬ 
ing H 2 {kr) by H 2 {kr) + 4:i/TT{kr)^, thus suppressing the 
near helds. The corresponding eigenvalue distribution is 
shown in Figj^b). Despite the fact that the ±1 sublevels 
remain coupled, the eigenvalue distribution of the vecto¬ 
rial case without near held terms closely resembles that 
of the scalar case, even though the smallest values of 7 „ 
do not reach the lowest limits obtained in the scalar case. 
Conversely, if we add a near held term, which we choose 
as the one present in vectorial scattering —4i/7r(fcr)^, to 
the scalar kernel HQ{kr), the long lived modes of the 
purely scalar case disappear [see FigH c)]. The scaling 
analysis, as well as a thorough analysis of the spatial ex¬ 
tension of the modes, conhrm that Anderson localization 
is absent from these altered interactions. 


Focusing on lifetimes, the study of the longest of 
them Tmin first reveals that for low densities {p < 0.3), 
long lifetimes are caused by the radiation trapping: 
~ l/6g [Hj, with 6o the cloud optical depth [15] (see 
Figj^c)). However, for scalar light, the appearance of the 
localized modes for p/k^ > 0.3 comes along with lifetimes 
much larger than those predicted by radiation trapping, 
see Figj^c). These lifetimes are not simply a function 
of the density but appear to scale as N'^/^/kR and 
to decay exponentially fast (see Figj^d)). This result is 
clearly beyond the standard Anderson localization, where 
quantities scale as N/R?, or cooperative effects where it 
scales as N/R, and calls for new approaches. Finally, 
while scalar light with near-fields exhibits lifetimes that 
always decay as 1 /6g (it is almost with the radiation trap¬ 
ping black curve in Figj^c)), both vectorial light with 
and without near-held exhibit lifetimes longer than that 
of radiation trapping: These come from atom pairs in¬ 
stead of localized modes, as reveals the analysis of the 
IPR and of the spatial profiles. 

These results suggest that both the presence of near 
field interaction terms [4| and coupling of different sub- 
levels can break down long lifetimes and localization. We 
also found that removing the ani sotro py present in vec¬ 
torial scattering in Eqs. 18b ) does not restore 

localized modes. 


Our 2D study, apart from the investigation of subra¬ 
diance and localization in lower dimensions, allows for 
a more efficient numerical study of the eigenvectors of 
the dipole-dipole coupling. One aspect of the eigenvec¬ 
tor analysis is already seen in Figj^ where the IPR of 
the eigenmodes allows, for instance, a clear identification 
of atomic pairs (red circles in Figcorresponding to an 
IPR close to 0.5, indicating atom pairs). In addition the 
2D configuration allows for an easy systematic study of 


the shapes of the eigenvectors: two typical eigenmodes 
are shown in Fig|^ The localized mode is spatially well 
confined [insert in FigQa)] and has a clear exponential 
shape over several orders of magnitude [Figj^a)]. Vec¬ 
torial eigenvectors, on the other hand, are extended over 
almost the whole system size [insert in Figj^b)], with no 
indication for an exponential decrease [FigQb)]. This 
observation is again in line with previous conclusions in 
3D |U [TU] . The scalar light with near fields and vectorial 
light without near fields do not exhibit any exponentially 
localized modes, but rather extended modes, as can be 
observed in Fig.[^ These scatterings, as well as the vec¬ 
torial light with near fields, may however present features 
of hybrid states where localized and extended subradiant 
features combine |I6j . 



10 '‘ lo” lo' 10 ^ 10 ^ 10 '^ 10 "' 10 '‘ , io“ io‘ 

p/k' 


FIG. 4. (Color online) Spatial profile of the most subradi¬ 
ant non-pair mode for (a) scalar light (exponentially localized 
mode) and (b) vectorial light (delocalized mode), as a func¬ 
tion of the distance to its center of mass Tcm. The 2D profile is 
exhibited in the inset, (c) Inverse lifetime versus localization 
length of the modes for scalar light; the localization length ^ is 
obtained by exponential fit of the spatial profile of the mode 
((a) for an example), so it is meaningful only for localized 
modes, within the dashed-bordered box. More specifically, 
the strongly subradiant modes lying outside of the box are 
extended, and so are the superradiant (^ typically exceeds 
the system size for these, sign of a full delocalization), (d) 
Localization length for scalar light as a function of the normal¬ 
ized density, for different scatterer number. The ’theory’ line 
refers to the theoretical prediction j4p) exp(7rfe^/8p). 

Panels (a-b) were realized for N = 5000 and p/k^ = 1, so 
kR ~ 40; Panel (c) is for N = 5000 and p/k^ = 10, so 
kR ~ 12.6, as marked by the dash-dotted line. 

The localized nature of the strongly subradiant modes 
is thus confirmed by the analysis of their spatial profile. 
Furthermore, as long as the mode does not considerably 
extend over the edge of the atomic cloud, its localization 
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FIG. 5. (Color online) Spatial profile of the most subradiant 
non-pair mode for (a) scalar light with near fields and (b) 
vectorial light without near helds. Simulations realized for 
an homogeneous disk cloud of = 5000 particles with an 
homogeneous p/k^ = 1 density 

length and lifetime are uncorrelated (uniform filling of 
the box in Figj^ with modes). There is however a corre¬ 
lation between the position of the mode (indicated by the 
color code, where blue points mark modes at the center 
of the system) and its localization length. Some modes 
even mix to surface (whispering gallery) modes that have 
a much larger spatial extend. Yet there is no correlation 
between the position of the mode (at the center or near 
the edge of the cloud) and its lifetime. The absence of 
correlation between the lifetime of the modes and their 
localization length calls for a differentiation between spa¬ 
tial and temporal localization. Although all spatially lo¬ 
calized modes are subradiant, the shortest localization 
length may not be associated to the longest lifetimes. 
This corroborates studies on photon escape rates that 
failed to observe the localization phase transition [7] , and 
is also highlighted by the fact that spatial localization is 
affected by boundary effects while temporal localization, 
surprisingly, does not seem to be. Finally, as can be 
seen in FigQd), for densities above ~ 0.05 the lo¬ 
calization length no longer depends on the system size, 
but only on the spatial density. These curves are not in 
agreement with the prediction of localization length from 


the perturbative approach in the weak disordered regime 

m- 

We have also verified (see Figj^a)) that the corre¬ 
sponding spectral overlap g and the scaling function f3 
for the altered interactions are qualitatively similar to 
the one of the purely vectorial case, i.e., they deviate only 
slightly from zero. Similarly, inspection of the eigenval¬ 
ues did not reveal any spatially localized mode. Together 
with the above results on lifetimes, this observation sug¬ 
gests that extremely long lifetimes of modes, well beyond 
radiation trapping ones, come along with spatial local¬ 
ization, i.e., subradiance may be a condition necessary to 
localization. 

IV. CONCLUSIONS 

In conclusion, we explored 2D scattering by point scat- 
terers in a scalar and a vectorial limit. Even though 
our eigenvalue analysis is consistent with previous results 
and interpretations of localization, our procedure of ar¬ 
tificially introducing or removing near field terms com¬ 
bined to a spatial analysis of the eigenfunctions support 
that very long lifetimes come along with Anderson lo¬ 
calization, but both near-held terms and the coupling of 
polarizations may prevent their emergence. Furthermore, 
we reported an absence of correlations between lifetime 
and localization length of localized modes, pointing at 
the difference between spatial and temporal localization. 
An important task for the future will be to relate both 
the 2D and 3D studies to transport properties of electro¬ 
magnetic waves and to compute observables that can be 
tested in experiments. 
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